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ABSTRACT 

We study the physics of electron acceleration at coUisionless shocks that move through a plasma 
containing large-scale magnetic fluctuations. We numerically integrate the trajectories of a large 
number of electrons, which are treated as test particles moving in the time dependent electric and 
magnetic fields determined from 2-D hybrid simulations (kinetic ions, fluid electron). The large- 
scale magnetic fluctuations effect the electrons in a number of ways and lead to efficient and rapid 
energization at the shock front. Since the electrons mainly follow along magnetic lines of force, the 
large-scale braiding of field lines in space allows the fast-moving electrons to cross the shock front 
several times, leading to efficient acceleration. Ripples in the shock front occuring at various scales 
will also contribute to the acceleration by mirroring the electrons. Our calculation shows that this 
process favors electron acceleration at perpendicular shocks. The current study is also helpful in 
understanding the injection problem for electron acceleration by coUisionless shocks. It is also shown 
that the spatial distribut ion of energetic electrons is similar to in-situ observations (e.g.. iBale et al.l 
[TQQilSimnett et al.l[2005h . The process may be important to our understanding of energetic electrons 
in planetary bow shocks and interplanetary shocks, and explaining herringbone structures seen in 
some type II solar radio bursts. 

Subject headings: acceleration of particles - cosmic rays - shock waves - turbulence 



1. INTRODUCTION 

CoUisionless shocks are widely believed to be the pri- 
mary acceleration mechanism giving rise to the ubiq- 
uitous existence of energetic particles in space. The 
theory of diffusive shock acceleration (DSA ) was pro- 
posed some 30 years ago (lAxford et al.lll978t lBellllI978l: 
IBlandford fc Ostrikedfl978t IKrvmskvl 119771) and is cur- 
rently believed to be the most important mechanism for 
a variety of astrophysical environments, for example, in- 
terplanetary shocks, the heliospheric termination shock, 
and shocks associated with supernova remnants. The 
theory predicts a universal power-law energy flux spec- 
trum dJ/dE oc for strong shocks with a density 
compression ratio of 4. 

The physical mechanism by which particles are accel- 
erated from thermal energies to much higher energies 
where DSA is presumed to be applicable (the injection 
problem) has received some recent attentions but has 
not reached a common consensus explanation. Many ac- 
celeration theories, for example, shoc k drift acceleration 
(SDA) (see reviews. [Armstrong et"aLiri9 85: Decker 1988} 



and shock surfing acce leration ( Sagdeevlil966t ILee et al.i 
[1996: Za nk et al lll996[ ) have been proposed. The accel- 
eration of low energy protons in the shocks containing 
large- scale pre-existing magnetic fluctuations is very effi- 
cient (|Giacalonell2005al [bl). which suggests that there may 
not be an injection problem. 

It is generally thought that, pre-accelerated particles 
will interact resonantly with magnetic turbulence which 
results in isotropization and diffusion. Many previ- 
ous works co nsidering magnetic turbulence focused on 
the ions (e.g.. lBelllll 978^ 'Giacal one et~al][l992[ : iNg et al.l 
[2003t IGiacalond 1200 4, 2 005a . b). However, the acceler- 
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ation of electrons is less well understood since for elec- 
trons whose gyroradii are very small, the cyclotron res- 
onance condition is not easily satisfied thus they can- 
not interact resonantly with large-scale ambient turbu- 
lence on ion-sc ale. While the scat tering provided by 
whistler waves (iShimada et al.l Il999f ) is one possibility, 
iJokipii fc Giacalond (j2007l ) proposed an attractive solu- 
tion to the injection problem that does not require pitch- 
angle scattering, i.e., conserving the first adiabatic invari- 
ant. The idea is that the low-rigidity particles, especially 
electrons, can move rapidly along meandering magnetic 
field lines and thus travel back and forth between shock 
front. The particles gain energy from the difference be- 
tween upstream and downstream flow velocities. 

Energetic electrons are often observed to be associ- 
ated with coUisionless shocks. Accelerated electrons 
are thought to produce type II radio b ursts in the so- 
lar co rona and interplanetary space. [Anderson et al.l 
(|1979f l reported ISEE spacecraft measurements of up- 
stream electrons (> 16 keV) of the Earth's bow shock 
that originate from a thin region close to the point of 
tangency between i nterplanetary magnetic fi eld lines and 
the shock surface. [Tsurutani fc LinI (|1985l ) showed ob- 
servations of energetic electrons associated with inter- 
planetary shocks showing " spike-like" flux enhancements 
for energies > 2 kcV. The spike events were observed at 
quasi-perpendicular shocks with 0Bn ^ 70°, where 0b„ 
is the angle between upstream magnetic field and shock 
normal. Some shock crossings had no enhancements of 
energetic electrons which were reported t o be associated 
with low shock speeds and small Osn- ^Si mnett et al.l 
(|2005[ ) presented data which shows energetic electrons 
are accelerated close to shock front. They also showed 
some accelerated electrons can escape far upstream of 
quasi- perpendicular interplanetary shocks. The clear ev- 



2 



Guo and Giacalone 



idence of electron acceleration b y PSA is rare, but a 
recent example was discussed by IShimada et al.l (|19990 
showing the importance of whistler wayes. 

In order to expl ain th e energ izatio n of electrons within 
the s hock layer, iWul ()1984[ ) and iLeroy fc Mangenevl 
(jl984[ ) developed analytic models for electron acceler- 
ation from thermal energies by adiabatic reflection by 
a quasi-perpendicular shock. This is known as fast- 
Fermi acceleration. This theory describes a scatter-free 
electron acceleration process in a planar, time-steady 
shock. It obtains a qualitative agreement with ob- 
servations at Earth's bow shock in terms of the loss- 
cone pitch-angle dist ribution and energy r ange of ac- 
celerated electrons. iKrauss-Varban et al.l ()1989D used 
the combination of electron test particle simulation and 
1-D hybrid simulation and verified Wu's basic conclu- 
sions. The main energy source of fast Fermi accelera- 
tion comes from t he —V x B/c electric field which is 
the same as SDA (jArmstrong et al1ll985[ ). It can also 
be demonstrated that fast-Fermi acceleration and SDA 
are t he same process in tw"0 diff erent frames of refer- 
ence (jKrauss-Varban fc Wul [19891 ). Thus one would ex- 
pect electrons to drift in the direction perpendicular to 
the fiow and magnetic field. For a single refiection, the 
fracti on and energies of acc ele rated particles are limited 
(e.g., IBall fc Melrosi [200l . iHolman fc PessesI ([l98l 
proposed the basic outline for type II solar radio bursts in 
which energetic electrons are accelerated through SDA. 
It is expected that multiple reflections are required in or- 
der to explain herringbone structures in type II bursts, 
where the electrons are acceler ated to a fractio n of the 
speed of light. More recently, iBurgesd ()2006[) studied 
electron acceleration in 2-D quasi-perpendicular shocks 
using test-particle simulations and self-consistent hybrid 
simulations. He focused on the effect of the rippling of 
the shock front on particle acceleration in highly oblique 
shocks with dsn > 80°. The ripples, in this case, were 
produced by instabilities along the shock front. Burgess 
found that the acceleration of electron by SDA can be 
more efficient for a rippled shock. The shape of the re- 
sulting energy spectra has a fiat plateau from the initial 
release energies to energies several times higher than this. 
Above the flat plateau the spectra drop off steeply as 6'b„ 
become smaller. 

In this paper, we use test particle simulations combined 
with 2-D hybrid simulations that in clude pr e-existing 
large-scale magnetic fleld turbulence (jGiacalon e 2005bi) 
to study the shock acceleration of electrons. In addi- 
tion to the effect of the large-scale turbulence, the shock 
microphysics occuring on ion length and time scales is 
also included. Moreover, the shock front is rippled and 
distorted in response to the turbulence, which is also in- 
cluded in our model. In § 2 we describe the numerical 
method we used to combine the fields from the hybrid 
simulation and test-particle simulation to obtain the elec- 
tron distribution. § 3 gives the main results of our sim- 
ulation. In § 4 we summarize the main conclusions and 
discuss the implication of our work. 

2. NUMERICAL METHOD 

Investigating particle transport in the vicinity of a col- 
lisionless shock requires a spatial scale large enough for 
particles to propagate back and forth across the shock, 
and a spatial resolution small enough to include the 



detailed physics for particle scattering and shock mi- 
crostructure. We implement a combination of a 2-D 
hybrid simulation to model the fields and plasma fiow 
and a test particle simulation to follow the orbits of a 
large number of energetic electrons. In the first step, 
we employ a two- dimensiona l hybrid simulation similar 
to previous work ([Giacalond l2005b[ ) that includes pre- 
existing la rge-scale turbulence. In the hybrid simula- 
tion (e.g., IWinske &: Questl Il988[ ). the ions are treated 
fully kinetically and thermal (i.e., non-energetic) elec- 
trons are treated as a massless fluid. This approach is 
well suited to resolve ion-scale plasma physics which is 
critical to describe supercritical coUisionless shocks. In 
this study, we consider a two-dimensional Cartesian grid 
in the x — z plane. All the physical vector quantities 
have components in three directions, but depend spa- 
tially only on these two variables. A shock is produced 
by using the so-called pis ton method (for a discussion, 
see iJones fc EllisonI Il991t ) , in which the plasma is in- 
jected continuously from one end {x = 0, in our case) 
of the simulation box, and reflected elastically at the 
other end (x = L^:). This boundary is also assumed 
to be a perfectly conducting barrier. The pileup of den- 
sity and magnetic fleld creates a shock propagating in 
the —x direction. To include the large-scale magnetic 
fluctuations, a random magnetic fleld is superposed on 
a mean field at the beginning of the simulation and is 
also injected continuously at the a; = boundary dur- 
ing the simulation. The simplified one-dimensional fiuc- 
tuations have the form B(z,t) = S'B{z,t) + Bi, where 
Bi is the averaged upstream magnetic field. The fluc- 
tuating component contains an equal mixture of right- 
and left-hand circularly polarized, forward and backward 
parallel-propagating plane Alfven waves. The amplitude 
of the fluctuations is determined from a Kolmogorov 
power spectrum: 



in which Lc is the coherence scale of the fluctuations, see 
(|Giacalonell2005bl) for more details. For the simulations 
presented in this study, we take Lc ~ L^, which is the 
size of simulation box in z direction. Note that in ad- 
dition to magnetic fluctuations, there are also velocity 
perturbations with 5v = vai6'B/Bi (Alfven waves). For 
most of the parts in the paper, we consider a turbulence 
variance a = 5B^/B\ = Sv^ /v\^ = 0.3, where Sv and 
VAi are the magnitude of velocity perturbation and up- 
stream Alfven speed, respectively. We also discuss the 
effect of different values of turbulence variances. 

The size of the simulation box L^- x Lz = iOOc/ojpi x 
1024c/wpi. where c/upi is the ion inertial length. The 
Mach number of the flow in the simulation frame is 
Mao = 4.0, the averaged Mach number in the shock 
frame is about 5.6. Most of the results presented here 
are for averaged shock normal angle < Osn >= 90°, 
but we also simulate the cases for < Osn >= 60° and 
75° to examine the dependence of the acceleration ef- 
ficiency on shock normal angle. The other important 
simulation parameters include electron and ion plasma 
beta ~ 0.5 and = 0.5, respectively, grid sizes 
Ax = Az = O.Sc/cjpi, time step At = 0.01fl~^, the 
ratio between light speed and upstream Alfven speed 
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c/vAi — 8696.0; and the anomalous resistivity is taken 
to be 77 = 1 X 10^^47rw~/. The initial spatially uni- 
form thermal ion distribution was generated using 40 
particles per cell. Different from previous studies, the 
consideration of large-scale magnetic fluctuations enables 
us consider the effect of pre-existing magnetic turbu- 
lence on electron acceleration, which has bee n shown to 
be impo rtant for low-energy ion acceleration (jGiacalond 
I2005aibr i since particle transport normal to the mean 
field is enhanced. However, the particle transport in full 
3-D turbulence can not be properly treated in a self- 
consistent way using availabl e computation. As demon- 
strated by previous works ([ Jokipii. Kota fc Gi acalond 
[19931 iGiacalone fc JokipiilllQQl Tjones et al. 1998), in the 
model with at least one ignorable coordinate, the center 
of gyration of particles is confined to within one gyrora- 
dius of the original magnetic field line. The test-electrons 
can still move normal to the mean field in our model be- 
cause of the field-line random walk. 

In the second part of our calculation we integrate the 
full motion equation of an ensemble of test-particle elec- 
trons in the electric and magnetic fields obtained in the 
hybrid simulations (see Figure 1). This part of the calcu- 
lation is done separately from the main hybrid simulation 
as a post processing phase. We assume non-relativistic 
motion which is reasonable because the highest energy 
electrons ob tained in our study are sti ll non-relativistic. 
As noted bv iKrauss-Varban et al.l pJ89ll . high-order in- 
terpolation of fields is required to ensure numerical accu- 
racy and avoid artificial scattering in calculating electron 
trajectories. In this work we use second-order spatial in- 
terpolation and linear temporal interpolation, which en- 
sure the smooth variations of the electromagnetic fields. 
We release a shell distribution of electrons with energy 
of 100 eV, which corresponds to an electron velocity 
Ve ~ 30.7L/i = b.lvthe in the upstream frame, where 
C/i is upstream bulk velocity in the shock frame and vthe 
is the thermal velocity of fluid electrons considered in the 
hybrid simulations, respectively. This energy is typical 
for the halo component of electron velocity distributions 
observed in solar wind. The test-particle electrons are re- 
leased uniformly upstream sX t = TOfij/ when the shock 
has fully formed and is far from the boundaries. The nu- 
merical technique used to integrate electron trajectories 
is the so-cal led Bulirsh- S toer m ethod, which is described 
in detail bv I Press et al.l (|1986[ ). It is highly accurate and 
conserves energy well. It is fast when fields are smooth 
compared with the electron gyroradius. The algorithm 
uses an adjustable time-step method based on the eval- 
uation of the local truncation error. The time step is 
allowed to vary between 5 x 10~^ and , where Vice 

is the electron gyrofrequency. The ratio Vlce/^d is taken 
to be the realistic value 1836. The total number of elec- 
trons in the simulation is 1.6 x 10^. The electrons which 
reach the left or right boundary are assumed to escape 
from the shock region and are removed from simulation. 
The boundary condition in the z direc tion is taken t o be 
periodic. The readers are referred to (|Burgesj|2006[ ) for 
more details on the numerical methods. 

Magnetic field turbulence has already been proved to 
have key effect on the particle acceleration in collision- 
less shocks. Unfortunately, solving the whole problem in 
three-dimensional space and resolving magnetic turbu- 



lence from coherence scale to electron scale are still lim- 
ited by available computation in the near future. This 
limitation motivates us to solve these problems approxi- 
mately. We also note that in our model the electron test 
particle simulation is not self-consistent since the hybrid 
simulation does not include the electron scale physics. 
The electron scale shock structure, which may be impor- 
tant is neglected here. 

3. SIMULATION RESULTS 
Figure 1 shows a snapshot of the z component of the 
magnetic field, B^/Bi, at time llOfi^/ in a gray-scale 
representation, overlayed by a two-dimensional magnetic 
vector field. At this time, the shock is fully developed. 
In this case, the angle between the average magnetic field 
direction and shock normal, < 0b„ > is 90°. The posi- 
tion of the shock front is clearly seen from the bound- 
ary of the magnetic field jump. The shock is moving 
in the —x direction at a speed dependent on z, which 
is about 1.6 vai on average. Because of the effect of 
large-scale turbulence with the shock, the shock surface 
become irregular on a variety of spatial scales from small- 
scale ripple s, which could be due to ion-scale plasma in- 
stabilities ()Lowe fc Burgess! 120031) , to large-scale struc- 
ture caused by the int eraction between the shock and 
upstream turbulence (iNeugebauer fc Giacalond 120051 : 
IGiacalone fc Neugebauedl2008t iLu et al Jl2Q0gD. The up- 
stream magnetic field is compressed and distorted as it 
passes through the shock into the downstream region. 
We note that the rippling of the shock and varying up- 
stream magnetic field leads to a varying local shock nor- 
mal angle along the shock front. As we will discuss later, 
the irregular shock surface and magnetic field geometry 
will efficiently accelerate electrons and produce a number 
of features similar to observations, such as the electron 
foreshock and spike-like intensity increases at the shock 
front. The meandering of field lines close to the shock 
surface helps to trap the electrons at the shock, leading to 
efficient acceleration. The shock ripples also contribute 
to the acceleration by mirroring electrons between them. 

Figure 2 shows a color-coded representation of the 
number of energetic electrons with energies higher than 
10 times (i.e., 1 keV) the initial (at release) energy at 
three different times (a) 76^-^, (b) 8117"/, (c) 90f2j/, 
respectively. It is found that after the initial release, a 
fraction of the electrons are reflected and accelerated at 
the shock front, and then travel upstream along the tur- 
bulent magnetic held lines. These accelerated electrons 
are then taken back to the shock by the field line me- 
andering, which provides even further acceleration. The 
number of energetic electrons close to the shock surface is 
highly irregular because the acceleration efficiency varies 
along the s hock front depending on the local shock nor- 
mal angle (|Wul I1984D . Most of the electrons are con- 
centrated near the shock front since the global magnetic 
field is mostly perpendicular to the shock normal. As 
the field lines convect through shock, the electrons even- 
tually are taken downstream. Since the electrons are 
tied to individual field lines in 2-D magnetic field, once 
the electrons are no longer capable of crossing the shock, 
there will be no additional significant acceleration. At 
this point, once all electrons are downstream, the energy 
spectrum no longer changes with time. 
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Examination of the trajectories of some electrons 
shows that the ripphng of the shock front also contributes 
to the acceleration by mirroring electrons between the 
ripples, as illustrated in Figure 3. In this figure, the top 
left plot displays the trajectory of a representative elec- 
tron in the x—z plane, overlapped with the 2-D gray-scale 
representation of at ftdt = 89.0, the gray scale is the 
same as in Figure 1. The upper right plot shows the po- 
sition of this electron (in x) as a function of time. The 
electron bounces back and forth between the ripples for 
several times. For example, the refiections are labeled 
a — b, c, d, f — g, h, and j. The energy change as a func- 
tion of position, x, corresponding to these reflections is 
shown in the bottom left panel. We find that there are 
jumps in energy at each of the reflections. The panel on 
the bottom right shows the electron energy as a function 
of time which also illustrates the features of multiple ac- 
celerations related to multiple reflections. The trajectory 
analysis shows the electron will be mirrored between the 
ripples for a couple of times and by this get accelerated 
multiple times. Note that the shock does not move much 
during the time scale of this trajectory. 

We now consider the effect of varying the angle be- 
tween the mean magnetic fleld and shock-normal. Shown 
in Figure 4 are the resulting energy spectra for three dif- 
ferent mean shock-normal angles (< Osn >— 60°, 75° 
and 90°, respectively) at the end of simulations {ftdt = 
120.0). It is found that for < 63,1 >= 90°, the electrons 
can readily be accelerated to up to 200 — 300 times the ini- 
tial energy within 50r2~^. The spectrum is flat between 
about 0.1 keV to 0.7 keV. This shape is similar to the 
"plateau" structure discussed bv lBurgessI ()2OO60 . Above 
1 keV, the spectrum falls off with energy as a slope index 
about —3. It can be found that both the number fraction 
and highest energy of accelerated particles decreases as 

< &Bn > decreases. We have also tried different value of 
initial energies (not shown), and find that the accelera- 
tion efficiency decreases for electrons wit h higher initia l 
energies, which is similar to the results of (lB urgessll2006D . 

The effect of different strengthes of magnetic turbu- 
lence is examined in Figure 5. We compare three cases 
with different turbulence variances a = 0.1, 0.3, and 0.5, 
respectively. At the end of simulations, the final energy 
spectra are similar at low energies, with significant vari- 
ations in the spectra only at energies larger than 2 keV. 
It is found that the energy spectrum is hardened at high 
energies when the turbulence variance is largest, which 
indicates the large-scale turbulence is more important for 
accelerating electrons to high energies. We argue that 
collisionless shocks which move through magnetic turbu- 
lence with significant power leads to efficient electron ac- 
celeration to high energies since the motion normal to the 
shock front is enhanced. The reason is that the meander- 
ing of field lines is enhanced, which allows the electrons 
have a better chance to travel though the shock multiple 
times. 

We note that the spatial distribution of energetic elec- 
trons is determined not only by the ripples in the shock 
front, but also by the global topology of the magnetic 
field lines. An example is shown in Figure 6, which 
shows the profiles of the number of energetic electrons 
at r^czi = 100.0 as a function of x, for the case of 

< (^Bn >= 90°. The black solid line is the profile at 



z = 200c/a;pi, and the red dash line shows the profile at 
z = SOOc/wpi. The corresponding position of the shock 
front at each of these values of z are represented using 
dot lines. At 2; = 200c/a;pi, it is observed that the ener- 
getic electrons travel far upstream up to about 100c/ LUpi. 
However, the profile at z = SOOc/tUpi shows no signifi- 
cant upstream energetic electron flux. The upstream en- 
ergetic electron profiles show irr egular features similar to 
in-situ observations reported bv lSimnett et al.l (|2005f ) (in 
Figure 10). The irregular features are controlled by the 
global topology of the large-scale turbulent magnetic field 
lines, along which the accelerated electrons could travel 
far upstream. Additionally, energetic electron profiles in 
X direction generally show " spike-like" structure close to 
the shock front, which is usually observed in interplane- 
tary shocks and Earth's bow shock. We note this feature 
is relatively stable within the simulation time once the 
upstream electron structure developed. 

4. DISCUSSION AND CONCLUSIONS 

We studied the acceleration of electrons at collisionless 
shocks by utilizing a combination of a 2-D hybrid simu- 
lation to obtain the shock structure and a test-particle 
simulation to determine the motion of electrons. The hy- 
brid simulation provides realistic electric and magnetic 
fields within the transition layer of the shock that effect 
the motion of test-electrons, which is determined by solv- 
ing the equation of motion. The interaction of the shock 
with pre-existing upstream fluctuations, and other non- 
linear processes occuring in the hybrid simulation lead 
to a "rippling" of shock surface which also effects the 
transport of the electrons. 

We find that the electrons are efficiently accelerated 
by a nearly perpendicular shock. The turbulent mag- 
netic field, leads to field-line meandering which allows 
the electrons to cross the shock front many times. The 
rippling of the shock front also contributes to the acceler- 
ation by mirroring electrons between the ripples. For the 
case that the averaged shock normal angle < Obu >= 90° 
and turbulence variance a = 0.3, the electrons can read- 
ily be accelerated to up to 200 — 300 times the initial 
energy. The resulting spectrum is flat between about 0.1 
keV to 0.7 keV. At higher energies, the spectrum falls 
off with energy like a power law with a spectral slope of 
about —3. This acceleration process is more efficient at 
perpendicular shocks. As < Obu > decreases from 90° , 
both the number fraction and highest achievable energy 
of accelerated particles decreases. Based on our calcu- 
lations, we conclude that perpendicular shocks are the 
most important for the acceleration of electrons. The 
current study is helpful in understanding the injection 
problem for electron acceleration by collisionless shocks. 
It is also found that different value of magnetic turbu- 
lence variances strongly affects the maximum energy at- 
tainable. The case with larger turbulence variance has 
a flatter energy spectrum than the case of smaller tur- 
bulence variance, which suggests the enhanced motion 
of electrons normal to the shock front, due to enhanced 
fleld-line random walk, is of importance for the acceler- 
ation of electrons to high energies. 

In addition, we also found that the energetic elec- 
tron density upstream and downstream of collisionless 
shocks show filamentary structures (Figure 2). This 
could help explain electron spike-like events observed up- 
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stream a nd downstream of terr e strial and interplanetary 
shocks (lAnderson et all 119791: iTsurutani fc Lid 119851: 
iSimnett et al.ll2005[ ). Observation by Voyager 1 at the 
termination shock and in the heliosheath also show the 
evidence of electron spike- like enhancements at the shock 
front (.Decker et al.. .2005[ ) . The upstream spatial dis- 
tribution of energetic electrons shows irregular features 
which depend on both the irregularity in the shock sur- 
face and the global topology of magnetic field lines. At 
first the electrons are accelerated and reflected at the 
shock front, and then they travel upstream along the 
magnetic field lines. The electrons could be taken far 
upstream by field line random walk. This result can pos- 
sibly lead to an interpretation to the complex electron 



foreshock events recent ly observed to be associated with 
interp l anetary shocks (pale et al.l 119991 iPulupa fc Bald 
[200l . iBaieetall (119991) and lPulupa fc Bald (|2008n pro- 
posed the complex upstream electron events are resulted 
from large-scale irregularities in shock surface. In this 
paper we have demonstrated that the upstream electron 
flux may be controlled by both an irregular shock surface 
and large-scale meandering magnetic field lines. 
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Fig. 1. — A snapshot of magnetic field in z direction Bz/Bi represented in gray-scale at time flOn^^^, where Bi is the averaged upstream 
magnetic field strength. A two-dimensional vector field is also overlapped, which indicates the direction of magnetic field. The shock surface 
is shown to be rippled and irregular in different scales. 
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Fig. 2. — The number of energetic electrons with energies E > 10-Eoi the initial release energy Eg = 100 eV, at (a) flcit = 76, (b) 
Qi-it = 81, and (c) Qcit = 90, respectively. Initially electrons are released uniformly upstream at Qcit = 70. 
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Fig. 3. — A typical electron trajectory analysis which shows acceleration by multiple mirroring between ripples. The top left panel 
displays the trajectory of the representative electron in x — 2: plane, overlapped with contour of Bz magnetic field where the gray-scale is 
the same as that in Figure 1; The top right panel shows the position of the electron in x coordinate as a function of time; The bottom 
left panel illustrates the energy of the representative electron E/Eq as a function of x; The bottom right panel shows the dependence of 
electron energy E/Eq on time. 
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Fig. 4. — The energy flux spectrum of electrons at Qcit = 120 for different averaged shock normal angle. The red solid line is in the case 
that the shock angle < 9g„ >= 90°, the blue dot dashed line and the black dashed line arc in the cases that < 9g„ >= 60° and 75°, 
respectively. 
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Fig. 5. — The energy flux spectrum of electrons at Qcit = 120 for averaged perpendicular shock with different turbulence variances. The 
black dashed line, red solid line, and blue dot dashed line are in the cases that a = 0.1, 0.3, and 0.5, respectively. 
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Fig. 6. — The solid black line and dashed red line show the profiles of the number of energetic electrons at 2 = 200c/ uipi and z = 800c /ujpi 
at time Qcit = 100, respectively. The red dot line and black dot line label the corresponding positions of the shock fronts. 



